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Abstract 

We point out that the measurements of the CP asymmetries in neutral B- 
decays together with a measurement of Br(K^ — > -r°vv) and the known value 
of | V us | can determine all elements of the Cabibbo-Kobayashi-Maskawa matrix 
without essentially any hadronic uncertainties. In particular the strong depen- 
dence of Br{K L — > ix°vv) on | V cb | (| V c b | 4 ) can be used to determine this element 
precisely. Consequently | V u b |, | Vtd | and | V ts \ can also be determined. We 
derive a set of formulae which facilitate this study and we present a numerical 
analysis of several future scenarios. Accurate determinations of all parameters 
are expected from the B-factories and the Main Injector era at Fermilab including 
KAMI. An impressive accuracy should be reached in the LHC era. An analysis 
using Br(K + — > 7i + uu) instead of Bt(Kl — > ^vv) shows that although the de- 
termination of | V c b | and | V t d \ this way is less accurate, still useful results for 
these elements can be obtained. 
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The Cabibbo-Kobayashi-Maskawa matrix [I], ^ which parametrizes the weak 
charged current interactions of quarks contains four parameters which must be 
determined by comparing the theoretical branching ratios with the experimental 
data. One of these parameters, the element V us , is known from K + — > 7i°e + i>, 
Kl — > n~e + v and semileptonic hyperon decays with a high precision: | V us \ = 
A = 0.2205 ±0.0018 || Recent critical discussions of this determination and of 
the related element | V u & | can be found in || . The present determinations of the 
remaining three parameters are subject to theoretical uncertainties resulting from 
our inability to perform precise non-perturbative calculations of various hadronic 
matrix elements of weak currents and local four-quark operators which enter the 
relevant formulae. 

It has been stressed in numerous papers 0, [J| that CP asymmetries in B° — > 
/, where / is a CP eigenstate can determine two parameters q and rj in the 
Wolfenstein parametrization H without essentially any hadronic uncertainties. A 
recent analysis of the related unitarity triangle can be found in [[J , where a high 
accuracy of this determination at future B-factories has been emphasized. Here 
we would like to point out that the measurement of the purely short distance CP 



violating decay Kl — > t\°vv pi together with the CP asymmetries in question, 



allows a clean and precise determination of the fourth parameter | | (or A). 
The argument is as follows. 



The last year calculations |p| , [12|j of next-to-leading QCD corrections to 
Br(Ki —>■ 7r°z/z/) reduced the theoretical uncertainty due to the choice of the 
renormalization scales present in the leading order expressions [fTS] below ±1%. 
Because the relevant hadronic matrix element of the weak current 57^(1— 7s)<i can 
be measured to better than 1% in the leading decay K + — > TT°e + i>, the resulting 
theoretical expression for BriK^ — > 7[ uu) is only a function of 77, V c b and m t = 
r?Tt(mt), the running top quark mass at the m t scale. The very weak dependence 
on the QCD scale Ajfs can be safely neglected. Now as shown below rj will be 



known with high precision from future B-factories [14], HERA-B [15], Tevatron 
with the Main Injector [|16[ and LHC B-physics experiments |F7|, [18| . Similarly 
Tevatron and LHC will offer precise determinations of m t . Consequently | V c b \ 
can be determined by measuring BriK^ — ■> n^vv). Since in addition BriK^ — > 
-k^vv) is a very sensitive function of | V c b \ (| V c b | 4 ), the latter element of the 
CKM matrix can be extracted from Kl — > -n^vv with a precision comparable or 
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even better than achieved directly in tree level semi-leptonic B-decays. To this 
end Br(Ki — > ^vv) must be known to within 5%. Since the present standard 
model estimate of BriK^ — > -n^vv) is in the ball park of (20 — 40) ■ 10~ 12 , the 
sensitivity of 0(1O- 12 ) aimed for at KAMI |TJ and KEK should allow to 
achieve this goal. 

At first sight it is probably surprising that we use a rare K-meson decay to 
determine | |. The natural place to do this are of course tree level B-decays. 
On the other hand using unitarity and the Wolfenstein parametrization with 
| V c b |= AX 2 it is clear that | V c b | gives the overall scale A of the top quark 
couplings Vtd and Vt s which are the only CKM couplings in Kl — > tt i>v. From 
this point of view it is very natural to measure the parameter A in a short distance 
process involving the top quark and using unitarity of the CKM matrix to find 
the value of | V c b |. Moreover this strategy in contrast to tree- level B-decays is 
free from hadronic uncertainties. On the other hand one should keep in mind 
that this method contains the uncertainty from the physics beyond the standard 
model which could contribute to short distance processes like —> 7t°z/za We 
will return to this below. 

The aim of this letter is to present this strategy for a clean determination of 
g, i] and | V^> | in explicit terms and to calculate the accuracy for these three 
parameters which we think should be achieved in the first decade of the next 
millennium. As byproducts we predict very accurate determinations of the top 
quark couplings | V t d | and | V ts | and of the ratio | V u b/V c b | which in our 
opinion are superior to any other determinations. We derive a set of formulae 
which should be useful in future studies of CP asymmetries. We also present an 
analysis in which Br(K + — > 7i + uu) instead of Br{Ki — > 7i uu) is used. We find 
that although the determination of | V c b | and | Vtd \ this way is less accurate, 
still useful results for these elements can be obtained. 

Our discussion of the Cabibbo-Kobayashi-Maskawa matrix will be based on 
the standard parametrization ||21|| , which can equivalently be rewritten in terms 
of the Wolfenstein parameters (A, A, g, rf) through the definitions 

Sl2 = A s 23 = AX 2 s 13 e~ iS = AX 3 (g - ir]) (1) 

Due to the resulting simplifications, the Wolfenstein parametrization M is par- 
ticularly useful when an expansion in A = \V US \ = 0.22 is performed. Including 
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next-to-leading terms in A || implies that the apex of the reduced unitarity tri- 
angle defined through 

Q + ifji 



Vu d V u \ 



is with an error of less than 0.1% given by 



X 2 



A 2 



(3) 



and not by (g, rj) as usually found in the literature |§. Working in the Wolfenstein 
parametrization such a treatment is required if we aim at a determination of the 
CKM parameters with a high precision. 

The CP asymmetries in 5° s -decays to CP eigenstates determine sm(2<f) i ) 
where fa = a, [3, 7 are the angles in the unitarity triangle defined by (0). 
Strategies involving several channels are sometimes necessary in order to remove 
hadronic uncertainties ||. sin(20j) can be expressed in terms of (g, fj) as follows 



sin(2a) 



sin(2/5) 



2r](r] 2 + g 2 - g) 
g 2 + 7/ 2 )((l - g) 2 + rf) 

27/(1 - Q) 



sin(27) 



1 - g) 2 + rf 
2gfj 



Next p| 



g 2 + rf 

Br (K L -> n°uu) = R L rf \ V cb | 4 X 2 (x t ) 



(4) 
(5) 
(6) 

(7) 



where x t = m 2 /M 2 



w ■ 



3a 2 Br(K+ -» n°e+v) r{K L ) 

KL = _ „ . ^ „ 7TTTT" = 0.354 • 10 



and 



x 



27r 2 sin 4 e T y r(K+) 
2 + x 3x — 6 



+ 



In a; 



^ = 0.985 



(9) 



1 — x (1 — a:" 12 

Here r] X is the NLO correction calculated in O, [T2|. With m 4 = rni(rrit) the 
QCD factor r/x is practically independent of mt. 

Using (0-^) we can express g, fj and | | in terms of 



a = sin(2o;) 6 = sin(2/3) c = sin(27) 



(10) 
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and Bt(Kl — > tt uu). There are several solutions. We give first only the solution 
which is favoured on the basis of what we already know about the CKM matrix. 
Using ((§-§) we express g in terms of fj and a, b, and c respectively: 



r] 2 + r]r^(a) , g = 1 — 7]r + (b) , g = r/r_(c) (11) 



where we have introduced 



r±{z) 



1 



i ± vT 



a, b, c 



:i2) 



For the pairs (a,b), (b,c) and (a,c), assuming 77 7^ 0, we can then determine 
fj with the result 

r_(a)+r + (6) 1 r_(a)+r_(c) 



l + r 2 + (b) r + (6)+r_(c)" l + r 2 (c) 
respectively. Using (|7|) and (W) we can next determine | V c b \ as follows: 



(13) 



V, 



cb 



X 7 



r]X(x t ) 



1/2 



Bj 



Br{K L n°vv) 



1.94-10 



-10 



(14) 



where 77 is to be found from (^) and ([13]). Note that the factor in front of A 2 gives 
the parameter A in the Wolfenstein parametrization. Finally using 



X(x t ) = 0.65 • x 



0.575 
t 



(15) 



which reproduces the function X(x t ) to an accuracy of better than 0.5% for 
150 GeV < m t < 190 GeV we find a useful formula 

1/4 

(16) 



cb 



39.1 • 10" 



/0.39 



V 



170 GeV 



0.575 



Br(K L 



3-10 



-11 



We note that the weak dependence of | Vcb \ o n Br(Ki — > 7i uu) allows to achieve 
high accuracy for this CKM element even when BriK^ — > -n^vv) is known within 
5 — 10% accuracy. 

Equations (|ll])-(0) together with ([3D are the main formulae which we will 
use to determine p, 77 and | V c b \- First however we would like to discuss the 
remaining solutions. 

There are four solutions for p coming from @ at fixed sin(2a) and fj. They 
are given by the first formula in (|ll|) with r±(a) and with ± in front of the 
square root. The solutions with + in front of the square root can be excluded by 
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imposing | V u b/V c b |< 0.10 in accordance with the data on semi-leptonic B-decays 
. Next there are two solutions for p coming from @ given by the second 
formula in (|TT|) with r±(b). The solution with r_(6) violates | V^/V^, |< 0.10 and 
consequently there is only one acceptable solution coming from (|^). Finally there 
are two solutions for p from given by the last formula in ( |TT| ) with r±(c). 



Retaining the allowed solutions for p we find then the generalization of ([13D 
with (r ± (a), r + (6)), (r + (6), r ± (c)) and (r±(a), r ± (c)), in an obvious notation, re- 
spectively. One can then check numerically by varying (a, b, c) in the full range 
that the unique solution with (r_(a), r+(5)), (r + (6), r_(c)) and (r_(a),r^(c)) as 
given in (|TTD and (|13|) is obtained if one requires 

?7 > 0.20, -0.20 < q < 0.25, 0.06 <| |< 0.10 (17) 

Vcb 

This range is favoured by a recent analysis of the unitarity triangle and the 
data on B-decays |22| . We will discuss only this solution in what follows. In the 
future when (a, b, c) will be measured one will have to use a similar strategy to 
select a unique solution by means of other measurements. Once this has been 
done, a precise determination of CKM parameters within this solution will be 
possible along the lines discussed here. We will illustrate this on examples below. 

We now turn to a numerical investigation of the formulae above. In (|]) and 
( |f4|) we use [^1[ 

A = 0.22 a = 1/128 sin 2 Q w = 0.23 (18) 

Br(K + -f 7r°e + z/) = 4.82 ■ 10" 2 t(K l )/t(K + ) = 4.18 (19) 

and assume that the 0(1%) uncertainties in these numerical values will be reduced 
in the coming years to the level that they can be neglected. We will also neglect 
the small residual scale ambiguity in X(xt) JR] which can effectively be taken 



into account by introducing an additional error Am t < ±1 GeV. 

As illustrative examples, let us consider the following three scenarios: 

Scenario I 

sin(2a) = 0.40 ± 0.08 sin(2/3) = 0.70 ± 0.06 

Br{K L ^7c°vv) = (3.0 ±0.3) • lO" 11 m t = (170 ± 5) GeV 

Scenario II 

sin(2a) = 0.40 ± 0.04 sin(2/3) = 0.70 ± 0.02 

Br{K L -»• ttVp) = (3.00 ± 0.15) • 10- 11 m t = (170 ± 3) GeV 
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Central 


I 


II 


III 


p 


0.072 


±0.040 


±0.016 


±0.008 


V 


0.389 


±0.044 


±0.016 


±0.008 


1 v cb | /10- 3 


39.2 


±3.9 


±1.7 


±1.3 


V u b/V c b 


0.087 


±0.010 


±0.003 


±0.002 


I V td | /10- 3 


8.7 


±0.9 


±0.4 


±0.3 


I ^ | /lO" 3 


38.4 


±3.8 


±1.7 


±1.3 


sin(27) 


0.35 


±0.15 


±0.07 


±0.04 


B K 


0.83 


±0.17 


±0.07 


±0.06 




200 


±19 


±8 


±6 



Table 1: Determinations of various parameters in scenarios I-III with Fb\JBb 
given in MeV. The errors have been symmetrized so that sometimes the central 
values do not exactly correspond to the central values of the input parameters. 



Scenario III 



sin(2a) 

-> 7T°Vi?) 



0.40 ±0.02 
(3.00 ±0.15) • 10~ n 



sin(2/3) 
m t 



0.70 ±0.01 



(22) 



Br{K L -+ = (3.00 ± 0.15) • lO" 11 m t = (170 ± 3) GeV 
The accuracy in the scenario I should be achieved at B-factories [|14| , HERA-B 



31| , at KAMI Jl9| and at KEK |20| . Scenarios II and III correspond to B-physics 



at Fermilab during the Main Injector era [|16| and at LHC [[H], |l8j respectively. 
At that time an improved measurement of Br(K^ — > 7i uu) should be aimed for. 
The values of mt assumed here are in the ball park of the most recent results of 



the CDF collaboration p3fl . Since in accordance with the QCD corrections in [IJ 
we use here the current top quark mass at the scale m t , our values correspond to 



m 



phys 



(177 ± 3) GeV. This should be compared with m 



phys 
t 



(174 ±16) GeV 



reported by CDF |23|. 

The results that would be obtained in these scenarios for p, rj, \ V c b |, 
I V u b/V c b |, | Vtd I) I Vts I an d sin(27) are collected in table 1. To this end we have 
used 

V ub 



V, 



id 



v cb i \y/(i -ey+fj 2 



v cb 



(23) 



and the standard expression for | V ts \ |2T|. 
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Central 


V 


IV 


nr 


p 


0.000 


±0.016 


±0.007 


±0.004 


V 


0.342 


±0.048 


±0.017 


±0.009 


1 v cb | /10- 3 


41.2 


±4.5 


±1.9 


±1.4 


V u b/V c b 


0.075 


±0.010 


±0.004 


±0.002 


I V td | /10- 3 


9.6 


±0.9 


±0.4 


±0.3 


I ^ | /lO" 3 


40.2 


±4.4 


±1.8 


±1.4 


sin(2a) 


0.60 


±0.12 


±0.05 


±0.03 


B K 


0.69 


±0.15 


±0.06 


±0.05 




174 


±15 


±6 


±5 



Table 2: Determinations of various parameters in scenarios P-III'. 

In table 2 we show the results for the corresponding scenarios P-III' in which 
sin(27) is used instead of sin(2a). Here we set 

f 0.60 ±0.06 (P) f 0±0.08 (P) 

sin(2/3) = I 0.60 ± 0.02 (IP) sin(2 7 ) = J ± 0.04 (IP) (24) 

[ 0.60 ±0.01 (IIP) [ 0±0.02 (IIP) 

Furthermore we use m t = 180 GeV and BriK^ —>■ ^vv) = 3.2 • 10" 11 with the 
errors as in scenarios I-III. We are aware of the fact that the accuracy for sin(27) 



assumed in scenarios IP and IIP will be difficult to achieve P, I24J but we show 



the results of this exercise to further motivate the efforts in this field. For p and 
7] we give three decimal places in order to show better the differences between 
various cases. One can easily check that other possible solutions for (g, if) can 
be excluded. Choosing for instance the solution (r + (a),r + (b)) in scenarios I-III 
one finds (g, r/) w (—1.6,1.1) and | V u b/V c b |~ 0.43 which is much too large. 
Similarly the solution (r+(6),r+(c)) in scenarios P-IIP gives (g,rj) ~ (1-0,0) and 
Vub/Vcb |~ 0.22 which is also too large. In the same manner all solutions 
discussed above except for the one presented in the tables can be excluded. 

Tables 1 and 2 show very clearly the potential of CP asymmetries in B-decays 
and of Kl — > -n^vv in the determination of CKM parameters. It should be stressed 
that this high accuracy is not only achieved because of our assumptions about 
future experimental errors in the scenarios considered but also because sin(2a) 
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and sin(27) are very sensitive functions of g and rj |J and Bt{Kl — > tt uu) 
depends strongly on | V c b |. In particular as seen in table 2 the determination of 
sin(27) may offer a very precise measurement of g. 

The accuracy in scenarios I and I' corresponds roughly to the cases consid- 
ered in ||, where however Kl — > tt°uu has not been discussed. There also the 
prospects of the determination of the unitarity triangle using B d — B d mixing and 
the parameter Ek have been analyzed. The results of show that in such an 
analysis it will be very difficult to determine g and rj to better than Ag = ±0.10 
and Arj = ±0.05 which should be contrasted with the accuracy expected here. 
Similarly because of theoretical uncertainties it is at present difficult to imagine 
that in the tree level B-decays a better accuracy than A | V c b \= ±2 • 10~ 3 and 
A | V u b/V c b |= ±0.01 could be achieved. 

In tables 1 and 2 we have also shown the values of the non-perturbative 
parameters Bk and Fb\/B~b which can be extracted from the data on Ek and 
B° d — B d mixing once the CKM parameters have been determined in the scenarios 
considered. To this end Xd = 0.72 and r(B) = 1.5 ps have been assumed. The 
errors on these two quantities should be negligible at the end of this millennium. 
Note that the resulting central values for Bk in tables 1 and 2 are close to the 



lattice |25| and l/N [f26j results respectively. 

We have checked that similar patterns of uncertainties emerge for different 
central input parameters. Needless to say when the scenarios presented here will 
become a reality one will have to make sure that the uncertainties present in the 
input parameters in fli"8]) and (fB5|) have been reduced to the desired level. 



It is instructive to investigate whether the use of another short distance decay 



K + — > n + uv instead of Kl — > it vis would also give interesting results for V c 



cb 



and Vtd- K + — > 7i + uu is CP conserving and receives also contributions from 
internal charm exchanges. This introduces the dependence on and m c but 
otherwise this decay is known to be theoretically very clean. In particular the long 



distance contributions to K + — > ix+vv have been considered in |27], g^, |29[ and 
found to be very small: two to three orders of magnitude smaller than the short 
distance contribution at the level of the branching ratio. Moreover in contrast 
to Kl — > 7T°z/z/, the decay K + — > 7i + uu could be observed already in the coming 
years at AGS in Brookhaven. Using the expressions in |30| and || we find instead 
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of (0) 



where 



V, 



cb 



1/2 

(77 2 + 6?)X(x t ) 



(25) 



and 



6! = 1 - g + A 2 (2f? - £ 2 - v 2 - i) + 0(A 4 ) (27) 

Here Pq(K + ) represents the charm contribution to K + — > n + uu calculated includ- 
ing next-to-leading QCD corrections in |3T||. For 200 MeV < < 350 MeV 
1.25 GeV <m c < 1.35 GeV, one has P (^ + ) = 0.40 ± 0.09 J3l], |3(J where also 
the residual uncertainty due to the choice of the renormalization scale /i has been 
taken into account. Here m c stands for the running charm quark mass at the m c 
scale. 

We again consider scenarios I-III with Br(K + — > 7r + z/z/) = (1.0 ± 0.1) • 10~ 10 
for the scenario I and Br(K + — > 7i + uu) = (1.0 ± 0.05) ■ 10~ 10 for scenarios II 
and III in place of Bt(Kl — > vr°z/z/) with all other input parameters unchanged. 
The results for g, r/, \ V u b/V c b \ and sin(27) remain of course unchanged. In 
table 3 we show the results for V c b, V t d and F B ^B B . We observe that due to 
the uncertainties present in the charm contribution to K + — > n + uu, which was 
absent in K L — ► tPvv, the determinations of | V^, |, | Vtd I an d F b ^Bb are 
less accurate. On the other hand, we note in the case of | V t d \ that due to 
an approximate cancellation of the square root in (^3p by the denominator in 
(|25"1), the uncertainties due to the errors in r\ and g are considerably reduced. 
Consequently when compared with — > -n^vv the reduction of the accuracy 
in V t d is less pronounced than in the case of V c b. If the uncertainties due to the 
charm mass and are removed one day, only the uncertainty related to fi 
will remain in P Q (K + ) giving AP (K + ) = ±0.03 pi[ . In this case the results in 
parentheses in table 3 would be found. 

Let us finally summarize the main aspects of this letter. 
We have pointed out that the measurements of the CP asymmetries in neutral 
B-decays together with a measurement of BriK^ — > -n°vv) and the known value 
of | V us | offer a precise determination of all elements of the Cabibbo-Kobayashi- 
Maskawa matrix without essentially any hadronic uncertainties. Kl —>■ ^vv 
proceeds almost entirely through direct CP violation and is known to be a very 
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Central 


/ 


II 


III 


v cb | /10- 3 


41.2 


±4.3 (3.2) 


±3.0 (1.9) 


±2.8 (1.8) 


v td /10- 3 


9.1 


±0.9 (0.7) 


±0.6 (0.4) 


±0.6 (0.4) 


Fb^/Bb 


190 


±17 (12) 


±12 (7) 


±12 (7) 



Table 3: Determinations of various parameters in scenarios I-III using K + 
7i + uu instead of Kl — > -n^vv. The values in parentheses show the situation when 
the uncertainties in m c and A-rje are not included. 



useful decay for the determination of rj. However due to the strong dependence on 
V c b this determination cannot fully compete with the one which can be achieved 
using CP asymmetries in B-decays. As a recent analysis shows it will be 
difficult to reach Arj = ±0.03 this way if | V c b | is determined in tree level B- 
decays. Our strategy then is to find r\ from CP asymmetries in B decays and use 
Kl — > -r^vv to determine | V c b |. To our knowledge no other decay can determine 
| V c b I as cleanly as this one. 

We believe that the strategy presented here is the theoretically cleanest way to 
establish the precise values of the CKM parameters. The ratio Xd/x s of B^ — B^ 
to B° s — B° s mixings, Br(K + — > 7r + z/p) and the parity violating asymmetry (Alr) 
in K + — > 7r + /i/2 |32|, are also theoretically rather clean and are useful 

in this respect. In particular as demonstrated in [[HJ, K + — > 7i + uu together 
with Kl — > -r^vv offers a respectable determination of sin(2/?). However from 
the theoretical point of view the determinations of this type cannot compete 
with the strategy considered here. Xd/x s is subject to uncertainties related to 
SU(3) flavour breaking effects which will probably be difficult to bring below 5%. 
Moreover the dependence on A cancells in this ratio and consequently V cb cannot 
be determined. Br(K + — ► 7i + uu) is more sensitive to Aj^g than Br(Ki — > r^vv) 
and in addition receives a few % uncertainty from the error in the charm quark 
mass pi |. As we have seen above these uncertainties lower the precision of the 
determination of V c b and Vu compared to Kl —>■ 7r°z/z/, although as shown in table 
3 interesting results for these elements can still be obtained. Similar comments 
can be made about A^^ for which in addition possible long distance contributions 
at a few % level cannot be excluded [|33l. 
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On the other hand once p, rj and | V c b | (or A) have been precisely determined 
as discussed here, it is clear that Xd/x s , Br(K + — > 7i + uu) and Alr can be 
rather accurately predicted and confronted with future experimental data. Such 
confrontations would offer excellent tests of the standard model and could possibly 
give signs of a new physics beyond it. 

Of particular interest will also be the comparision of | V c b | determined as 
suggested here with the value of this CKM element extracted from tree level 
semi-leptonic B-decays. Since in contrast to Kl — > 7t ujj, the tree-level decays 
are to an excellent approximation insensitive to any new physics contributions 
from very high energy scales, the comparison of these two determinations of | Vcb \ 
would be an excellent test of the standard model and of a possible physics beyond 
it. Also the values of | V u b/V c b | from tree-level B-decays, which are subject to 
hadronic uncertainies larger than in the case of V c b, when compared with the 
clean determinations suggested here could teach us about the realiability of non- 
perturbative methods. The same applies to the quantities like Xd and the CP 
violating parameter Ek which are subject to uncertainties present in the non- 
perturbative parameters Fb\/Bb and Bk respectively. In tables 1-3 we have 
shown the values of these parameters which can be extracted from the data in 
the scenarios considered. 

It is also clear that once the accuracy for CKM parameters presented here has 
been attained, also detailed tests of proposed schemes for quark matrices |}5, |36| 
will be possible. 

Precise determinations of all CKM parameters without hadronic uncertain- 
ties along the lines suggested here can only be realized if the measurements of 
CP asymmetries in B-decays and the measurements of Bt(Kl — > is^vv) and 
Br(K + — > 7i + uu) can reach the desired accuracy. All efforts should be made to 
achieve this goal. 

I would like to thank Gerhard Buchalla, Robert Fleischer and Markus Laut- 
enbacher for critical reading of the manuscript and Markus Lautenbacher for 
checking my numerical analysis. I also thank Hubert Kroha and Karl Jakobs for 
illuminating discussions about future B physics experiments. 
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